Structural, mechanical, thermodynamic, and electronic properties of thorium hydrides 

from first principles 
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We perform first-principles calculations of the structural, electronic, mechanical, and thermody- 
namic properties of thorium hydrides (TI1H2 and TI14H15) based on the density functional theory 
with generalized gradient approximation. The equilibrium geometries, the total and partial densities 
of states, charge density, elastic constants, elastic moduli, Poisson's ratio, and phonon dispersion 
curves for these materials are systematically investigated and analyzed in comparison with exper- 
iments and previous calculations. These results show that our calculated equilibrium structural 
parameters are well consistent with experiments. The Th— H bonds in all thorium hydrides exhibit 
weak covalent character, but the ionic properties for TI1H2 and TI14H15 are different due to their 
different hydrogen concentration. It is found that while in TI1H2 about 1.5 electrons transfer from 
each Th atom to H, in TI14H15 the charge transfer from each Th atom is around 2.1 electrons. 
Our calculated phonon spectrum for the stable body-centered tetragonal phase of TI1H2 accords 
well with experiments. In addition we show that ThH2 in the fluorite phase is mechanically and 
dynamically unstable. 
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I. INTRODUCTION 

Thorium is one kind of important nuclear materials 
and together with its compounds has been widely investi- 
gated both experimentally and theoretically. Among tho- 
rium compounds, thorium hydrides (TI1H2 and TI14H15) 
are metallic solids and have potential use for advanced 
nuclear fuels. In addition, TI14H15 has been reported to 
have superconductivity with transition temperatures T c 
^8 K 1 and has been considered as promising candidates 
for hydrogen storage since its large hydrogen-to-metal 
ratio 2 -. ThH 2 is not superconducting above 1 K£. 

Despite the abundant research on thorium hydrides, 
relatively little is known regarding their chemical bond- 
ing, mechanical properties, and phonon dispersion. Only 
the optical phonon density of states of TI1H2 and TI14H15 
was measured through inelastic neutron scatterings and 
the bulk modulus of TI1H2 was calculated by linear 
muffin-tin orbital (LMTO) method^. Until now the elas- 
tic properties, which relate to various fundamental solid- 
state properties such as interatomic potentials, equation 
of state, phonon spectra, and thermodynamical proper- 
ties, are unknown for ThH2 and TI14H15. In addition, 
although the electronic properties as well as the chemi- 
cal bonding in thorium hydrides have been calculated re- 
cently by Shein et al£ through the full-potential LAPW 
(FLAPW) method, the study of the bonding nature of 
Th— H bond involving its mixed ionic/covalent charac- 
ter is still lacking. These facts, as a consequence, in- 
hibit deep understanding of thorium hydrides. Motivated 
by these observations, in this paper, we present a first- 
principles study by calculating the structural, electronic, 
mechanical, and thermodynamic properties of thorium 



hydrides. Also, the stability of the metastable phase of 
TI1H2 (fluorite structure with space group Fm3m) is an- 
alyzed and our calculated results show that the fluorite- 
type ThH2 is mechanically unstable. We perform the 
Bader analysis^ of thorium hydrides and And that about 
1.5 (2.1) electrons transfer from each Th atom to H for 
ThH 2 (Th 4 H 15 ). 



II. COMPUTATIONAL METHODS 

Our total energy calculations are carried out by em- 
ploying the plane-wave basis pseudopotential method 
as implemented in Vienna ab initio simulation package 
(VASP)£. The exchange and correlation effects are de- 
scribed by the density functional theory (DFT) within 
generalized gradient approximation (GGA) 9 . The pro- 
jected augmented wave (PAW) method of BlochliS is 
implemented in VASP with the frozen-core approxima- 
tion. The thorium 6s 2 7s 2 6p 6 6d 1 5f 1 and the hydrogen 
Is 1 electrons are treated as valence electrons. Note that 
although the 5/ states are empty in elemental Th, this 
level turns to evolve into a hybridization with the hydro- 
gen orbitals both in the valence band and the conduc- 
tion band, as well as to prominently contribute to the 
conduction band (see Fig. [3] below) . 9x9x9 and 5x5x5 
Monkhorst-Packii k point-meshes in the full wedge of 
the Brillouin zone are used for ThH2 and Tt^His, re- 
spectively. Electron wave function is expanded in plane 
waves up to a cutoff energy of 450 eV, and all atoms are 
fully relaxed until the Hellmann-Feynman forces become 
less than 0.02 eV/A. 

In present work, the theoretical equilibrium volume, 
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bulk modulus B, and pressure derivative of the bulk mod- 
ulus B 1 are obtained by fitting the energy-volume data 
in the third-order Birch-Murnaghan equation of state 
(EOS) 12 . In order to calculate elastic constants, a small 
strain is applied onto the structure. For small strain e, 
Hooke's law is valid and the crystal energy E(V, e) can 
be expanded as a Taylor series^ 3 -, 



E(V, e) = E(V , 0) + V J2 + 



Vr 



(1) 



i,j=l 



where E(Vo,0) is the energy of the unstrained system 
with the equilibrium volume Vo, e is strain tensor which 
has matrix elements (i,j=l, 2, and 3) defined by 



ei 
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(2) 



and Cij are the elastic constants. For cubic structures, 
there are three independent elastic constants (G\i, C12, 
and C44). So, the elastic constants for fee TI1H2 and 
bec TI14H15 can be calculated from three different strains 
listed in the following: 



(5,5,5,0,0,0),e 2 = (5,0,(5,0,0,0), 
e 3 = (0,0,0,5,(5,(5). 



(3) 



The strain amplitude 5 is varied in steps of 0.006 from 
0.036 to 0.036 and the total energies E(V, 5) at these 
strain steps are calculated. After obtaining elastic con- 
stants, we can calculate bulk and shear moduli from the 
Voigt-Reuss-Hill (VRH) approximations^—. The Voigt 
(Reuss) bounds on the bulk modulus By (Br) and shear 
modulus Gy (Gji) for these two cubic crystal systems are 
deduced from the formulas of elastic moduli in Ref.— . As 
for bet TI1H2, the six independent elastic constants (Cn, 
C12, C44, C13, C33, and Cqq) can be obtained from six 
different strains listed in the following: 



(5,0,5,0,0,0), 
(0,0,5,0,0,0), 
= (0, 0, 0, 5, 5, 0), e 6 = (0, 0, 0, 0, 0, 5) 



(5, 5, 5,0,0,0),e 2 
(5,5,0,0,0,0),e 4 



(4) 



and the formulas of clastic moduli in VRH 
approximations— are from Ref. 15 . Based on Hill 
approximation^, B=^(Br + By) and G=^(Gr + Gy). 
The Young's modulus E and Poisson's ratio v are given 
by the following formulas: 

E = 9BG/{W + G),v= (3B - 2G)/[2(3B + G)]. (5) 




(b) 



FIG. 1: (Color online) Crystal structures of (a) bet TI1H2, 
(b) fee TI1H2, and (c) bee TI14H15. Here, larger cyan spheres 
stand for Th atoms and the smaller white H. Note that the 
magenta and red spheres in (c) stand for HI and the white 
and grey spheres H2. 
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FIG. 2: Comparison of total energy vs the primitive cell vol- 
ume for ThH2 in bet and fee phases. 



III. RESULTS 

A. Atomic and electronic structures of thorium 
hydrides 

At ambient condition, the stable thorium dihydridc 
crystallizes in a body-centered tetragonal (bet) ionic 
structure with space group IA/mmm (No. 139). Its unit 
cell is composed of two TI1H2 formula units with the tho- 
rium atoms and the hydrogen atoms in 2 a (in Wyck- 
off notation) and Ad sites, respectively [see Fig. 1(a)]. 
Each Th atom is surrounded by eight H atoms forming 
a tetragonal and each H connects with four Th atoms 
to build a tetrahedron. The present optimized lattice 
parameters (a and c) are 4.0670 A and 4.9125 A (see 
Table U), in good agreement with the experimental^ val- 
ues of 4.10 A and 5.03 A. Besides, the face-centered cu- 
bic (fee) fluorite-type structure with space group Fm3m 
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TABLE I: Optimized structural parameters, partial and total densities of states at the Feimi level N(Ef) (states/eV/Th atom) 
for thorium hydrides. For comparison, other theoretical results and available experimental data are also listed. 

Compound Cell constants (A) Coordinates Th 6d Th 5/ H s Total 

bet ThHa a=4.0670 (4.0475)°, (4.10) c Th (2a): 0, 0, 0.383 0.165 0.002 0.996 

(I4/mmm) c=4.9125 (4.9778)°, (5.03) c H (4d): 0, 0.5, 0.25 (0.207)° (0.140)° (0.001)° (0.866)° 

fee ThH 2 a=5.4851 (5.4902)°, (5.489)° Th (4a): 0, 0, 0.610 0.244 0.003 1.557 

(Fm3m) H (8c): 0.25, 0.25, 0.25 (0.374)° (0.202)° (0.002)° (1.451)° 

Th 4 Hi 5 a=9.1304 (9.1280)°, (9.11) c ' d Th (16c): x, x, x; z=0.2087 (0.208) d 0.610 0.678 0.031 2.334 

(I43d) HI (12a): 0.375, 0, 0.25 (0.441)° (0.619)° (0.049)° (2.474)° 

H2 (48e): 0.372, 0.219, 0.404 (0.4, 0.23, 0.372) d 

° Reference^, b Reference^, c Reference^, d Reference^. 
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FIG. 3: (Color online) Total and orbital-resolved local densities of states for (a) bet TI1H2, (b) fee TI1H2, and (c) bec TluHis. 
The Fermi energy level is set at zero. 



(No. 225) is considered as metastable phase for TI1H2 
[see Fig. 1(b)]. The fluorite structure is the stable struc- 
ture of all actinide dioxides. Moreover, this fee structure 
can be obtained from the bet stable structure by modu- 
lating the base vectors. In fee structure, each Th atom is 
surrounded by eight H atoms forming a cube and each H 
connects with four Th atoms to build a tetrahedron. Our 
optimized lattice constant (a) for fee ThH 2 is 5.4851 A 
(see Table H|, in excellent agreement with the measured^ 
values of 5.489 A. The structure of TI14H15 is more com- 
plicated. It crystallizes in a body-centered cubic (bec) 
structure with space group /43c! (No. 220) [see Fig. 1(c)]. 
The unit cell has sixteen Th (16c), twelve HI (12a), and 
forty eight H2 (48 e) atoms. Each HI is surrounded by 
four Th atoms to build the tetrahedral structure and each 
H2 is located at the center of triangle formed by three Th 



atoms. As a result, each Th is surrounded by three HI 
atoms and nine H2 atoms. In Fig. 1(c), we label clearly 
the HI atoms (red spheres) and H2 atoms (grey spheres) 
connected to the Th atom located at center. The current 
optimized structural parameters of TI14H15 are listed in 
Table HI where one can find that our calculation results 
are in good agreement with experiment^. 

To investigate the stability of ThH 2 , we have optimized 
the structural parameters of its bet and fee structures at 
different pressures. To avoid the Pulay stress problem, 
we perform the structure relaxation calculations at fixed 
volumes rather than constant pressures. For fee struc- 
ture, due to its high symmetry, the structure relaxation 
calculations are performed at fixed volumes with no re- 
laxation of coordinates and cell shape. However, for bet 
structure, the cell shape is necessary to be optimized due 
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to their internal degrees of freedom. The total energies 
(per primitive cell) of the two structures at different vol- 
umes are calculated and shown in Fig. [2] Obviously, 
the bet TI1H2 is more stable than fee TI1H2 under am- 
bient pressure. The equilibrium volumes of bet and fee 
structures are 40.63 A 3 and 41.26 A 3 , respectively. Thus 
the equilibrium volumes of both bet and fee phases are 
approximately equal and the total energy difference at 
their equilibrium states is ^0.032 eV. This conclusion is 
consistent well with the previous FLAPW results 2 -. 

Basically, all the macroscopical properties of materi- 
als, such as hardness, elasticity, and conductivity, origi- 
nate from their electronic structure properties as well as 
chemical bonding nature. Therefore, it is necessary to 
perform the electronic structure analysis of thorium hy- 
drides. The calculated total densities of states (DOS) and 
the orbital-resolved partial densities of states (PDOS) of 
bet Thfi2, fee ThET2 and bee TI14H15 are shown in Fig. 
[3l Overall, the occupation properties of ThH 2 in bet 
and fee phases are similar. The occupied DOS near the 
Fermi level is featured by the three well-resolved peaks. 
The one near —1.0 eV is principally Th 6c? in character, 
while the other two peaks respectively near —4.5 and 
—6.0 eV are mostly H Is hybridized with Th 6c? and 
Th 5/ orbitals. These well-separated orbital peaks have 
been observed in the photoelectron spectroscopy (PES) 
measurement 2 -. In addition, the band width of the Th 
6c? valence band near the Fermi level is 3 eV, consistent 
with the FLAPW calculation^ and experimental data 21 . 
The H Is valence band width is of 5.0 eV, also in ac- 
cord with the experimental data (6.0 eV) and previous 
FLAPW result (5.2 eV). Moreover, the low bands of bet 
(fee) ThH 2 covering from -18.8 (-19.0) to -17.2 (-17.5) 
eV is mainly featured by Th 6p state mixed with a little 
H Is state and the conduction bands is principally occu- 
pied by Th 5/ states with admixtures of the Th 60? and 
H Is states and has a width of 5.5 (5.1) eV. For TI14H15, 
our calculations reproduce all the features calculated by 
Shein et al£ and agree well with PES measurements 2 -. 
Note that our calculated conduction band width is 1.7 
eV, consistent with the PES experimental valued (0.9 
eV) and evidently smaller than that of FLAPW calcula- 
tion result* 2 -. The main occupation of conduction bands 
is Th 5/ orbital, mixed with a little Th 6c? and H Is 
states. In addition, we have also presented in Table|T]thc 
partial and total DOSs at the Feimi level N(Ep) for tho- 
rium hydrides to study the occupation of the conduction 
band. Clearly, our results are in good agreement with 
the FLAPW calculation 2 .. 

In order to gain more insight into the bonding nature 
of ground state thorium hydrides, we also investigate the 
valence charge density distribution. The calculated va- 
lence charge density maps of the bet ThH2 in (100) plane, 
fee TI1H2 in (110) plane, and bee TI14H15 in plane estab- 
lished by three Th and one H2 atoms are plotted in Fig. 
HI Clearly, the charge densities around Th and H ions are 
all near spherical distribution with slightly deformed to- 
wards the direction to their nearest neighboring atoms. 




(a) _ (b) " (c) 



FIG. 4: Valence charge density of (a) bet TI1H2 in (100) plane, 
(b) fee ThFb in (110) plane, and (c) bec TluHis in plane 
established by three Th and one H2 atoms. Contour lines are 
drawn from 0.0 to 1.0 at 0.05 e/A 3 intervals. 
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FIG. 5: (Color online) Line charge density distribution be- 
tween Th atom and the nearest neighbor H atom for bet 
ThH 2 , fee ThH 2 , and bec Th 4 Hi 5 . 



To describe the ionic/covalent character quantitatively 
and more clearly, we plot in Fig. [5]the line charge density 
distribution along the nearest Th— H bonds and perform 
the Bader analysis^. For bet and fee TI1H2, the two 
charge density curves change with the same way along 
the Th— H bonds. For TI14H15, the two curves along the 
Th— HI and Th— H2 bonds are almost the same before 
the bridge locus (indicated by the arrow), after where 
they split evidently. Besides, one can find a minimum 
value of charge density in Fig. [5] for each bond at around 
the bridge. These minimum values are listed in Table [TTJ 
Although these values are much smaller than 0.7 e/A 3 
for Si covalent bond, they are prominently higher than 
0.05 e/A 3 for Na— CI bond in typical ionic crystal NaCl. 
Also we notice that these values are smaller than the 
minimum values (0.45 e/A 3 ) of charge density along the 
Th— O bond in our previous study of ThO? 2 ^. There- 
fore, there are weak but clear covalent bonds between 
Th and H for TI1H2 in bet and fee phases and for bee 
Th 4 H 15 . The slight difference appears in Th 4 H 15 , where 
the minimum value of charge density in Th— HI bond is 
smaller than that in Th— H2 bond. In addition, one can 
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TABLE II: Calculated charges and volumes according to Bader partitioning as well as the Th— H distances and correlated 
minimum values of charge densities along the Th— H bonds for thorium hydrides. As a comparison, other theoretical work and 
available experimental data are also listed. Note that some properties of TI14H15 have two values (first for HI and second for 
H2). 



Compound 


Qs(Th) 


Qs(H) 


V s (Th) 


V S (H) 


Th-H 


Charge density m i n 




(e) 


(e) 


(A 3 ) 


(A 3 ) 


(A) 


(e/A 3 ) 


bet ThHa 


10.530 


1.735 


25.59 


7.50 


2.375 (2.385)" 


0.273 


fee ThEb 


10.488 


1.756 


25.72 


7.77 


2.375 (2.337)'' 


0.271 


Th 4 Hi 5 


9.870 


1.590, 1.563 


20.15 


7.16, 7.35 


2.465 (2.46)", 2.343 (2.29)" 


0.231, 0.296 



a Reference^, b Reference^. 



notice from Table [TT1 that our calculated Th— H distances 
are all well consistent with the previous theoretical^ and 
experimental^ results. As for Bader analysis, it is a 
well established analysis tool for studying the topology 
of the electron density and thus is particularly suitable 
for discussing the mixed ionic/covalent character of a 
compound. The charge (Qb) enclosed within the Bader 
volume (Vb) is a good approximation to the total elec- 
tronic charge of an atom. In the present study, the de- 
fault charge density grids for one unit cell are 56x56x72, 
56x56x56, and 112x112x112 for bet ThH 2 , fee ThH 2 , 
and bec Th4H 15 , respectively. To check the precision, 
the charge density distribution for fee ThH 2 is calculated 
with a series of n times finer grids (n — 2, 3, 4, 5, 6). The 
deviation of the effective charge between the five and the 
six times finer grids is less than 0.02%. Thus we perform 
the charge density calculations using the six times finer 
grid (336x336x432, 336x336x336, and 672x672x672 
for bet TI1H2, fee ThH 2 , and bee TI14H15, respectively). 
The calculated results arc presented in Table HU Note 
that although we have included the core charge in charge 
density calculations, since we do not expect variations as 
far as the trends are concerned, only the valence charge 
are listed. From Table [TT] the following prominent fea- 
tures can be seen: (i) The Bader charges and volumes for 
bet and fee ThH 2 are almost equal to each other. This 
shows the same ionic character, through a flux of charge 
(about 1.5 electrons for each Th atom) from cations to- 
wards anions, for thorium dihydride both in their sta- 
ble phase and metastable phase, (ii) The Bader charges 
and volumes for bec TI14H15 are different comparing to 
thorium dihydride. This originates from their different 
hydrogen concentration. For bec TI14H15, about 2.1 elec- 
trons transfer from each Th atom to H. This indicates 
that the Th atoms in TI14H15 are more ionized than that 
in thorium dihydrides. 



B. Mechanical properties of thorium hydrides 

Our calculated elastic constants, various moduli, pres- 
sure derivative of the bulk modulus B , and Poisson's 
ratio v for bet ThH 2 , fee ThH 2 , and bec Th4Hi5 are col- 
lected in Table IIIII Obviously, bet ThH 2 is mechanically 



stable due to the fact that its elastic constants satisfy 
the following mechanical stability criteria 2 ^ of tetragonal 
structure: 

Cn > 0, c 33 > 0, c 44 > 0, c 66 > 0, 

(C11 - C 12 ) > 0, (C n + C 33 - 2C 13 ) > 0, 

[2(Cn + C 12 ) + C 33 + 4C 13 ] >0. (6) 

TIi4Hi5 is also mechanically stable because its elas- 
tic constants satisfy the following mechanical stability 
criteria 2 ^ of cubic structure: 

C n > 0,C 44 > 0,C n > |Cia|,(C u +2Ci 2 ) > 0. (7) 

However, the ThH 2 in its metastable cubic phase is me- 
chanically unstable. In fact, one can see from Table 
IIIII that Cn is much smaller than Ci 2 for fee ThH 2 . 
Therefore, the mechanical stability criteria 2 ^ Cn>|Ci 2 | 
can not be satisfied and the shear modulus, Young's 
modulus, and Poisson's ratio also can not be obtained. 
As for bulk modulus B, each derived value from VRH 
approximations^ for all the thorium hydrides turns out 
to be very close to that obtained by fitting Murnaghan 
EOS 12 . Moreover, we notice that our calculated B for 
fee ThH 2 is approximately equal to previous LMTO the- 
oretical result^. For shear modulus G and Young's mod- 
ulus E, their values for bec TI14H15 are approximately 
2.5 times of those for bet ThH 2 . For Poisson's ratio, the 
values for both bet ThH 2 and bec TI14H15 are within the 
range from 0.25 to 0.45 for typical metals. 

C. Phonon dispersion curves of thorium dihydrides 

Employing the Hellmann-Feynman theorem and the 
direct method, we have calculated the phonon curves 
along some high-symmetry directions in the Brillouin 
zone (BZ), together with the phonon PDOS. For the 
phonon dispersion calculation, we use the 2x2x2 bet 
(fee) supercell containing 48 (96) atoms for bet (fee) 
ThH 2 . To calculate the Hellmann-Feynman forces, 
we displace four and eight atoms, respectively, for 
bet and fee ThH 2 from their equilibrium positions 
and the amplitude of all the displacements is 0.03 
A. The calculated phonon dispersion curves along the 
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TABLE III: Calculated elastic constants, various moduli, and Poisson's ratio v for thorium hydrides. 



Compound 


Cn 


Cl2 


C44 


C13 


C33 


Cee 


B 


B 


G 


E 


V 




(GPa) 


(GPa) 


(GPa) 


(GPa) 


(GPa) 


(GPa) 


(GPa) 


(GPa) 


(GPa) 


(GPa) 




bet ThH 2 


131 


100 


29 


78 


89 


5 


91 


3.3 


16 


44 


0.419 


fee ThH 2 


46 


119 


57 








95 (99)° 


3.8 








Th 4 Hi 5 


121 


66 


59 








85 


2.2 


43 


111 


0.282 



a Reference^. 




FIG. 6: (Color online) Calculated phonon dispersion curves and corresponding PDOS for (a) bet TI1H2 and (b) fee TI1H2. For 
comparison, the experimental data in Ref. [4] are also shown (squares). 



T-N-P-X-T-Z directions for bet ThH 2 and along 
the T-X-K-T-L-X-W-L directions for fee ThH 2 
are displayed in Fig. [BJa) and Fig. HJb), respectively. 

For both bet and fee TI1H2 , there are only three atoms 
in their primitive cells. Therefore, nine phonon modes 
exist in the dispersion relations. Due to the fact that 
thorium is much heavier than hydrogen atom, the vibra- 
tion frequency of thorium atom is apparently lower than 
that of hydrogen atom. Therefore, evident gap between 
the optic modes and the acoustic branches exits and the 
phonon DOS of bet (fee) ThH2 can be viewed as two 
parts. One is the part lower than 4.4 (4.1) THz where 
the main contribution comes from the thorium sublattice, 
while the other part range from 26.0 to 32.1 (26.4 to 31.8) 
THz is dominated by the dynamics of the light hydrogen 
atoms. In experimental measurements, Dietrich et al. 
reported that the acoustic branches of TI1H2 are lower 
than 4.8 THz and optic modes range from 23.4 to 36.2 
THa^. Their measured data of optic modes are presented 
in Fig. |5{a) for comparison and the acoustic-phonon data 
are unable obtained from their time-of-flight spectrum. 
Obviously, our calculated results are on the whole con- 
sistent with the experimental data. In addition, while 
all frequencies are positive for bet ThH 2 , which assures a 
dynamical stability of the bet phase against mechanical 
perturbations. For fee ThH 2 , one can see from Fig. E](b) 
that the transverse acoustic (TA) mode close to T point 
becomes imaginary along the T—K and T—L directions. 
The minimum of the TA branch occurs along the T—K 
direction. This indicates instability of fee phase of ThH 2 



compared to stable bet phase, which is well consistent 
with our previous mechanical stability analysis of ThH 2 . 
Moreover, there exists LO-TO splitting at T point in bet 
phase. But in fee phase there is no LO-TO splitting. 

IV. CONCLUSION 

In summary, we have used the first-principles DFT- 
GGA method to calculate the structural, electronic, me- 
chanical, and thermodynamic properties of ThH 2 in its 
stable bet phase and metastable fee phase and TI14H15 in 
bee phase. Our optimized structural parameters are well 
consistent with experiments. The occupation characters 
of electronic orbitals also accord well with experiments 
and previous calculations. Through Bader analysis, we 
have found that the Th— H bonds in all thorium hydrides 
exhibit weak covalent character, but the ionic property 
for ThH 2 and TI14H15 are different. While ~1.5 electrons 
transfer from each Th atom to H in ThH 2 , in ThjHis 
about 2.1 electrons deviate from each Th atom. In ad- 
dition, our calculated phonon curves of fee ThH 2 have 
shown that the TA mode becomes imaginary close to T 
point and the subsequent instability of this phase. 

Acknowledgments 

This work was supported by the Foundations for De- 
velopment of Science and Technology of China Academy 



7 



of Engineering Physics under Grants No. 2009B0301037 and No. 2009A0102005. 



* Author to whom correspondence should be addressed. E- 
mail: zhang_ping@iapcm.ac.cn 

1 C. B. Satterthwaite and I. L. Toepke, Phys. Rev. Lett. 25, 
741 (1970). 

2 I. R. Shein, K. I. Shein, N. I. Medvedeva, and A. L. 
Ivanovskii, Physica B 389, 296 (2007). 

3 C. B. Satterthwaite and D. T. Peterson, J. Less-Common 
Met. 26, 361 (1972). 

4 M. Dietrich, W. Reichardt, and H. Rietschel, Solid State 
Commun. 21 603 (1977). 

5 M. S. S. Brooks and B. Johansson, Physica B 130, 516 
(1985). 

6 R. F. W. Bader, Atoms in Molecules: A Quantum Theory 
(Oxford University Press, New York, 1990). 

7 W. Tang, E. Sanville, and G. Henkelman, J. Phys.: Con- 
dens. Matter 21, 084204 (2009). 

8 G. Kresse and J. Furthmiiller, Phys. Rev. B 54, 11169 
(1996). 

9 J. P. Perdew, K. Burke, and Y. Wang, Phys. Rev. B 54, 
16533 (1996). 

P. E. Blochl, Phys. Rev. B 50, 17953 (1994). 

1 H. J. Monkhorst and J. D. Pack, Phys. Rev. B 13, 5188 



(1972). 

12 F. Birch, Phys. Rev. 71, 809 (1947). 

13 J.F. Nye, Physical Properties of Crystals, Their Represen- 
tation by Tensors and Matrices, Oxford Press, Chap. VIII, 
1957. 

14 W. Voigt, Lehrburch der Kristallphysik (Teubner, Leipzig, 
1928). 

15 A. Reuss and Z. Angew, Math. Mech. 9, 49 (1929). 

16 R. Hill, Phys. Phys. Soc. London 65, 349 (1952). 

17 J. Hanies, J. M. Leger, and G. Bocquillon, Annu. Rev. 
Mater. Res. 31, 1 (2001). 

18 J. P. Watt and L. Peselnick, J. Appl. Phys. 51, 1525 (1980). 

19 C. Keller, Thorium (Springer, 1978). 

20 W. H. Zachariasen, Acta Cryst. 6, 395 (1953). 

21 J. H. Weaver, J. A. Knapp, D. E. Eastman, D. T. Peterson, 
and C. B. Satterthwaite, Phys. Rev. Lett. 39, 639 (1977). 

22 B. Wang, H, Shi, W. Li, and P. Zhang, aiXiv:0908.3558vl 
[cond-mat.mtrl-sci] (2009). 

3 J. F. Nye, Physical Properties of Crystals (Oxford Univer- 
sity Press, Oxford, 1985). 



